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Abstract 

Time-dependent density functional theory is thought to work well for the test cases 
of He and Be atoms. We perform a quantum defect analysis of the s, p, and d 
Rydberg states of Be with accurate ground state Kohn-Sham potentials. The s 
and p-quantum defects are well described by the ALDA, but fails badly for the 
d-quantum defect. The same failure is observed in case of He. This provides a new 
challenge for functional development in time-dependent density functional theory. 

1 Introduction 

Time-dependent density functional theory (TDDFT) [1] has become a popular 
method for calculating excitation energies of atoms, molecules, clusters, and 
solids [2,3]. In many common situations, it yields transition frequencies with 
errors below 0.2 eV, and properties of excited states are obtained with an 
accuracy comparable to that of ground-state DFT. 

As any theory is increasingly applied, the limitations of its approximations are 
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discovered [3]. An early limitation for TDDFT was due to the poor asymp- 
totic behavior of potentials in approximations to ground state DFT. This 
problem was overcome by the development of asymptotically corrected poten- 
tials [4,5,6,7]. Later, it was found that adding an asymptotic tail to the LDA 
potential, without any modification in the inner region, gives very accurate re- 
sults [8] . Another limitation was that within the adiabatic approximation, dou- 
ble excitations are missed by TDDFT. It was shown by Maitra et al. that this 
problem can be overcome with a frequency dependent exchange-correlation 
(xc) kernel [9,10]. A still open question is how to obtain charge-transfer exci- 
tations within TDDFT. An empirical approximation of the xc-kernel has been 
designed as a step to overcome this problem [11]. 

In this communication, we report a new and surprising failure of present 
TDDFT calculations. While previously hinted at in the hterature, our recently- 
developed quantum defect analysis (QDA) for atoms [12] makes it very clear 
that this is a qualitative failure. Essentially, as the angular-momentum change 
grows, TDDFT with the adiabatic local density approximation (ALDA), be- 
comes worse and worse. For s ^ d transitions, including TDDFT xc-contributions 
is much worse than ignoring them altogether! 



2 Theory 

For closed shell atoms and for any spherical one-electron potential that decays 
as — 1/r at large distances, the bound-state transitions form a Rydberg series 
with frequencies: 

1 

Unl^I - 7u (1) 
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where / is the ionization potential, and /ini is called the quantum defect [13]. 
We use atomic units {e^ = h = nie = 1) throughout. For real atoms, quantum 
defects depend only weakly on the principal quantum number n for large n 
and converge to a finite value in the limit n ^ oo. 

The quantum defect is a smooth function of energy, and is well approximated 
a polynomial of some low order p: 

li^\E)=j2fiiE\ E^u-I. (2) 

1=0 

Thus we calculate transition frequencies, extract the quantum defect via Eq. 1 
and plot them as a function of E. When we fit our quantum defect values to 
Eq. 2, we optimize the fit over the entire range of excitation energies, not just 
around E — Q (the //j are simply related to the a, 6, and c coefficients in 
Refs. [14] and [15]). 

For the TDDFT part, wc solve the linear response equations as derived by 
Casida [16]. Such a TDDFT calculation consists of two parts. First the ground 
state KS orbitals and KS orbital energies are determined with DFT. For He 
and Be, nearly exact potentials have been developed by Umrigar et al. (He [17] 
and Be [18]), and we use these potentials for the ground state. In this way, we 
eliminate any error due to the underlying ground state KS potential. In the 
second part of the calculation (the TDDFT part), we do perturbation theory 
on the ground state results. This means solving Casida's equations, and we 
use the ALDA xc-kernel in all cases. We also include results in the Hartree 
approximation, equivalent to setting the xc-kernel to zero. In this approxima- 
tion triplet values correspond to the bare KS orbital energy differences, as the 
Hartree interaction involves no spin fiips. 
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3 Computational details 



For the calculation of KS orbital energies and TDDFT excitation energies, we 



used the Amsterdam Density Functional program package ADF ( http: / /www.scm.com 
The numerical integration accuracy was set to 12 significant digits, and the 
convergence criterion in the self-consistent procedure for the solution of the 
KS equations was set to 10"^^. We reached convergence with a large even 
tempered basis set with primitive basis functions including / functions and a 
fit set including g functions. 

In order to map the accurate potentials onto the integration grid of ADF we 
performed a cubic spline interpolation between the available data points. We 
confirmed the quality of our ground state KS results by comparing the ADF 
orbital energies with orbital energies obtained with an OEP (optimized poten- 
tial model) program [19,20,21]. This program is basis set independent, works 
with a radial grid, and both the energies and the potentials are optimized in 
a self-consistent way. For this code, the accurate potentials do not need to be 
interpolated. We obtained ground state orbitals energies with ADF that were 
identical to the ones obtained with the OEP code at least up to 5 digits (in 
a.u.). 

In order to obtain the coefficients in Eq. 2 for He, we included states up to 
n = 6 in case of the ALDA and /xc = calculations, up to n = 10 in case 
of the accurate wave function reference data, and up to n = 7 for the KS 
values. For Be we included states up to n = 6 in case of ALDA and /xc = 
calculations, up to n = 8 in case of the experimental reference data, and up 
to n = 9 for the KS values. 
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Fig. 1. Hartree and ALDA singlet and triplet s-quantum defects of Be, compared 
with experimental values [23]. The closed symbols correspond to singlet values, the 
open symbols correspond to triplet values. 

4 Results 



Van Gisbergen et al. first calculated transitions from the accurate Be KS 
potential [22] using the ALDA kernel and reported tables with transitions. If 
the same data is cast into a quantum defect and plotted as function of energy, 
systematic errors are revealed and more insight in the results is obtained. We 
also did calculations on higher s — > d transitions that have not been studied 
previously. 

The continuous lines in all graphs we show correspond to the fitted quantum 
defect values. The coefficients used to obtain the continuous quantum defect, 
using Eq. 2, are reported in Ref. [12] for the s- and p-cases. We show the s- 
quantum defect of Be in Figure 1. The dashed line with open triangles shows 
the scattering from the ground state KS potential, which lies neatly between 
the experimental singlet (closed squares) and triplet (open squares) values. We 
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Be p-quantum defect 
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Fig. 2. Hartree and ALDA singlet and triplet p-quantum defects of Be, compared 
with experimental values [23]. The closed symbols correspond to singlet values, the 
open symbols correspond to triplet values. 

note how large the quantum defects are. We see from Eq. 1 that a quantum 
defect of 1 means (for single-particle scattering) that the deviation from — 1/r 
is so big that the energy of a particular transition to a state with quantum 
number n corresponds exactly to the energy of that transition in a hydrogen 
atom to a final state with quantum number n' — n — 1. Inclusion of only 
Hartree effects shifts this line down to the closed triangles for the singlet, but 
leave the triplets unchanged. Finally, the ALDA xc-kernel shifts both of these 
Hartree results upward to the circles. The ALDA accounts for more that 80% 
of the error of the Hartree approximation. 



We next show the the p-quantum defect in Figure 2. In this case we are seeing 
a lot more structure than in case of the s-quantum defect. We see that the 
ALDA slightly overestimates the curvature of the quantum defect correction 
in this case. This produces the largest errors for the most negative energies 
(i.e. the lowest transitions). When E goes to zero all quantum defect values get 
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Be d-quantum defect 
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Fig. 3. Hartree and ALDA singlet and triplet d-quantum defects of Be, compared 
with experimental values [23]. The closed symbols correspond to singlet values, the 
open symbols correspond to triplet values. 

close together except for the singlet Hartree values which underestimate the 
experimental results. The triplet Hartree results are a good approximation to 
the experiment. Here ALDA clearly overestimates the correction to /x(0). This 
is not a large problem for lower transitions, but causes the £■ — > results to be 
worse than bare Hartree. In fact, the KS values lie in-between the experimental 
singlet and triplet values until n = 6. It can be seen from the graph that they 
do seem to get below the experimental singlet and triplet results for n > 6. 
Indeed, for n = 7, both the singlet and triplet experimental quantum defects 
(0.3616 and 0.3676 respectively) lie above the KS value of 0.3550. So here 
we have an example where the quantum defect of the bare KS orbital energy 
differences does not lie between the experimental singlet and triplet values. 

We show the d-quantum defect of Be in Figure 3 and that of He in Figure 4. 
For Be we first note that the KS curve hes, as before, between the singlet and 
triplet. But in this case it is much closer, and parallel to, the triplet curve. In 
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He d-quantum defect 
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Fig. 4. Hartree and ALDA singlet and triplet d-quantum defects of He, compared 
with accurate wave function results [24] (labeled "Exact"). The closed symbols 
correspond to singlet values, the open symbols correspond to triplet values. 

this case the Hartree approximation does rather well, by shifting the singlet 
curve about halfway towards the experiment, while not moving the triplet. 
But the ALDA curves are a disaster. For the singlet, the ALDA correction is 
in the wrong direction! For the triplet the correction is in the right direction, 
but it is too large by a factor of about four. We find similar results in case 
of He (Figure 4). First of all, the KS curve lies below the singlet and triplet 
curves for all n values studied. Again it is close to the triplet curve, the singlet 
Hartree curve underestimates the reference curve. The ALDA correctly shifts 
both singlet and triplet curves upward, but the shift is too large for the singlet, 
again leading to a dramatic overestimation. 



We show the coefficients of the d-quantum defect corresponding to the fit of 
Eq. 2 in Table 1. The columns labeled KS and Hartree are both referring to 
the Hartree approximation, we label the triplet data by KS to emphasize that 
these values correspond to the bare KS orbital energy differences. Apart from 
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Table 1 



We give the coefHcients of the expansion of the singlet and triplet d-quantum defect 
of He and Be. The ionization energy of He is 0.9037 a.u., and that for Be is 0.3426 
in all cases. 







Reference ^ 


f 

J xc 


= 


ALDA 


Atom 




Singlet 


Triplet 


Singlet 


Triplet 


Singlet 


Triplet 


He 




0.0020 


0.0026 


-0.0034 


0.0013 


0.0349 


0.0087 










-0.0496 




u. iuoy 






Max. AE 


0.0002 


0.0005 


0.0006 


0.0003 


0.0002 


0.0003 


Be 




-0.0947 


0.1094 


0.0138 


0.0994 


0.1725 


0.1736 




Ml 


0.5636 


-0.0810 


0.2353 


0.0554 


0.5291 


0.1699 






-29.2439 


-1.1540 


7.8553 




20.0356 


1.7686 




Max. AE 


0.001 


0.001 


0.0003 


0.0002 


0.0000 


0.0008 



^ Accurate wave function results for He from Ref. [24] and experimental data for 
Be from Ref. [23]. 



the coefficients we also give the "Max. AE" in this tables, by which we mean 
the maximum absolute error between the fitted values and the true quantum 
defect values. In all cases we stopped adding coefficients to the fit expansion 
until this error was smaller than 0.001, or when adding more coefficients no 
longer reduced the error. We see that three coefficients are enough in all cases 
to completely describe the quantum defect. The ALDA estimate for the E = 
quantum defect (/io) is close to the exact data for Z = 0, 1, but not for / = 2 
and we also see that the E — quantum defect is not always between the 
exact singlet and triplet values. 
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Be d-quantum defect 
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Fig. 5. Hartree and TDOEP singlet and triplet d-quantum defects of Be, compared 
with experimental values [23]. The closed symbols correspond to singlet values, the 
open symbols correspond to triplet values. 

Thus in case of the d-quantum defect, including the ALDA kernel within 
TDDFT worsens the results instead of improving it. In the mid 1980s a similar 
failure of the LDA was observed in the case of transfer energies [25,26]. It was 
found that sp transfer energies in first row atoms and sd transfer energies in 
2>d atoms are underestimated by the local spin-density approximation (LSD). 
The source of the error was linked to the fact that the LSD is unable to 
properly take into account the nodal structure of the orbitals. The failure of 
the ALDA that we observe here may be related to the LSD problem, but the 
methodology here is very different. In TDDFT, we extract transitions by using 
linear response, while in order to obtain transfer energies, one subtracts total 
energy differences. 

To test the effect of including orbital dependence in the kernel we extracted 
quantum defect values from the exchange-only OEP calculations of Petersilka 
et al. [27]. They did not do full time-dependent OEP, instead they use the 
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approximate PGG kernel [28]. We show results in Figure 5. The errors are 
similar to those of ALDA, but considerably quantitatively smaller. A full OEP 
calculation might do better still. 

This work was supported by NSF Grant No. CHE-0355405. 
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